Anderson localization versus charge-density-wave formation in disordered electron systems 
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We study the interplay of disorder and interaction effects including bosonic degrees of freedom in the frame- 
work of a generic one-dimensional transport model, the Anderson-Edwards model. Using the density-matrix 
renormalization group technique, we extract the localization length and the renormalization of the Tomonaga 
Luttinger liquid parameter from the charge-structure factor by a elaborate sample-average finite-size scaling 
procedure. The properties of the Anderson localized state can be described in terms of scaling relations of the 
metallic phase without disorder. We analyze how disorder competes with the charge-density-wave correlations 
triggered by the bosons and give evidence that strong disorder will destroy the charge-ordered state. 

PACS numbers: 71.23.An, 71.27.+a, 71.30.+h, 71.45.Lr 



Disorder is an inherent part of any solid state system (]]]. 
Low-dimensional materials are exceedingly susceptible to dis- 
order. In one dimension (ID), theory predicts that all carri- 
ers are strongly localized for arbitrary energies and arbitrar- 
ily weak disorder. This holds for Anderson's non-interacting 
tight-binding Hamiltonian with a diagonal (i.e., on-site) ran- 
dom potential 12j,|3[]. The coherent backscattering from the 
randomly distributed impurities thereby transforms the metal 
into an insulator. 

In ID the mutual interaction of the particles is likewise of 
significance; here even weak interactions can cause strong 
correlations. The instantaneous Coulomb repulsion between 
the electrons, for instance, tends to immobilize the charge car- 
riers as well. As a consequence, at half-filling, a Mott insu- 
lating (spin-density-wave) phase is energetically favored over 
the metallic state |4]. The retarded electron-phonon coupling, 
on the other hand, may lead to structural distortions accompa- 
nied by polaron formation and is the driving force behind 
the metal-to-Peierls transition, establishing a charge-density- 
wave (CDW) order H. 

To understand how disorder and interaction act together is 
of vital importance not only discussing the metal-insulator it- 
self but also analyzing the electronic properties of many quasi- 
1D materials of current interest, such as conjugated poly- 
mers, organic charge transfer salts, ferroelectric perovskites, 
halogen-bridged transition metal complexes, TMT[SF,TF] 
chains, Qn(TCNQ)2 compounds, or, e.g., the quite recently 
studied vanadium dioxide Peierls-Mott insulator |7|4 1011 . Car- 
bon nanotubes ll 111 and organic semiconductors 11211 are other 
examples where disorder and bosonic degrees of freedom are 
of importance. Regarding interacting bosons, ultracold atoms 
trapped in optical lattices offer the unique possibility to tune 
both the disorder and interaction strength 1 1311 . 

Unfortunately the subtle interplay of disorder and interac- 
tion effects is one of the most challenging problems in solid 
state theory and — despite fifty years of intense research — still 
an area of uncertainty, see Il4ll and references therein. In 
the limit of vanishing charge carrier density only the inter- 
action with the lattice vibrations matters. Then Anderson dis- 
order may affect the polaron self-trapping in a highly nontriv- 



ial way II 1 511 . This has been demonstrated for the Anderson- 
Holstein model within the statistical dynamical mean field and 
momentum average approximations jig , Il7ll . At finite carrier 
density, the Mott- Anderson transition for Coulomb correlated 
electrons was investigated by self-consistent mean-field the- 



ory in D = oo and D = 3 II 1 8142011 . as well as by variational 
Gutz wilier ansatz based approaches 12111 . Electron-electron 
interactions may screen the disorder potential in strongly cor- 
related systems, stabilizing thereby metallicity ll22Tl . Exact re- 
sults are rare however. In 2D, Lanczos and quantum Monte 
Carlo data suggest a disorder-induced stabilization of the 
pseudogap, also away from half-filli ng 12311 . The density ma- 
trix renormalization group (DMRG) 12411 allows the numerical 
exact calculation of ground-state properties of disordered, in- 
teracting fermion systems in ID, on fairly large systems. Ex- 
ploiting this technique, the properties of disordered Luttinger 
liquids have been analyzed in the framework of the spinless 
fermion Anderson-i-V model (AtVM) 1E5I1 and the spinful 
Anderson-Hubbard model 12611. 



In this paper, we address how many-body Anderson local- 
ization competes with CDW formation triggered by bosonic 
degrees of freedom in the framework of the Anderson- 
Edwards model (AEM). The Edwards model \2% represents 
a very general two-channel fermion-boson Hamiltonian, de- 
scribing quantum transport in a background medium. Its 
fermion-boson interaction part 



H fb = -t b J2f]fM+ b j) 

(id) 



(1) 



mimics the correlations/fluctuations inherent to a spinful 
fermion many-particle system by a boson-affected transfer of 
spinless charge carriers. In Eq. (fl}, a fermion fy' creates 
(or absorbs) a local boson every time it hops to a nearest 
neighbor (NN) site j. Thereby it creates a local excitation in 
the background with energy ojo- Hb = <*>o J2i blK- Because 
of quantum fluctuations the background distortions should be 
able to relax with a certain rate A. The entire Edwards Hamil- 
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FIG. 1. (Color online) DMRG metal-insulator phase boundary for 
the ID half-filled Edwards model without disorder (solid line). CDW 
order is suppressed if the background fluctuations dominate [wq < 1] 
or if the system's ability for relaxation is high [A > X c (ojq)]. The 
blue crosses denote the parameter sets considered in this paper for 
the disordered Edwards model. 



tonian then reads 



H E =H fb -\J2(bl 



(2) 



A unitary transformation 6, i-> bi + X/u>o eliminates the boson 
relaxation term in favor of a second fermion hopping channel: 
H E = H fb - t f J2(ij) fj.fi + H b . We like to emphasize 
that (i) this free-fermion transfer takes place on a strongly 
reduced energy scale tf = 2Aib/wo however, and (ii) co- 
herent propagation of a fermion is possible even in the limit 
A = tf = by means of a six-step vacuum-restoring hopping 

The 



process |28|], acting as a direct next NN transfer "fj +2 fi"- 
Edwards model reveals a surprisingly rich physics. Depend- 
ing on the relative strengths tf/t b of the two transport mech- 
anisms and the rate of bosonic fluctuations tb/u>o it repro- 
duces Holstein and t-J model-like lattice- and spin-polaron 



transport, respectively, in the single -particle sector 1128 , 
For the half-filled band case, a metal-insulator quantum-phase 
transition from a repulsive Tomonaga luttinger liquid (TLL) 
to a CDW has been reported S QM, se e Fig- CD Note that 
the CDW is a few-boson state that typifies rather a corre- 
lated (Mott-Hubbard-type) insulator than a Peierls state with 
many bosons (phonons) involved (30[ 31]. Since in the limit 
uo > 1 > ) (here, and in what follows t b is taken as 
the unit of energy) background fluctuations are energetically 
costly, charge transport is hindered and an effective Hamil- 
tonian with NN fermion repulsion results. To leading or- 
der, in a reduced (zero-boson) Hilbert space, we get H t y = 



-*/E< 



///< 



t b /u)Q. This so- 



/ 

called t-V model can be mapped onto the exactly solvable 
XXZ model, which exhibits a TLL-CDW quantum phase 
transition at V/tf = 2, i.e., at A^ 1 = 4. This value is 

smaller than those obtained for the Edwards model in the limit 

l 







< 1, where A" 1 ~ 6.3 (see Fig. Q] and Ref. |3l||), be 



cause already three-site and effective next-NN hopping terms 
were neglected in the derivation of the tV-model. 

We now employ the DMRG technique 12411 . which can be 



easily generalized to treat systems including bosons 13211 . in 
order to obtain unbiased results for the full AEM, 



H 



AE 



e,n 



f 



(3) 



and the related AtVM, Haw = ^2~2i £ i n { + H w, where 
disorder of strength A is induced by independently distributed 
random on-site potentials e,-, drawn from the box distribution 
P(si) = 0(1/2 — \£i\). Within the pseudosite approach, a 
boson is mapped to n b pseudosites 11321 13311 . In the numer- 
ical study of the AEM we take into account up to n b — 4 
pseudosites and determine n b by the requirement that local 
boson density of the last pseudosite is less than 10~ 7 for all i. 
Furthermore we keep up to m = 1200 density-matrix eigen- 
states in the renormalization steps to ensure that the discarded 
weight is smaller than 10~ 8 . The calculations are performed 
for finite systems with lengths L = 16 to 128 and open bound- 
ary conditions (OBC). For the simpler effective AfVM we 
reach L = 192 with OBC. Here the use of m = 1000 density- 
matrix eigenstates makes the discarded weight negligible. To 
gain representative results for our disordered systems we pro- 
ceed as follows. We first compute the physical quantity of 
interest at fixed L for numerous samples {£;}, then set up an 
appropriate statistical average, and finally perform a careful 
finite-size scaling. 

An important question is, of course, which physical quan- 
tity to use in the finite-size scaling of the Anderson transi- 
tion. For this purpose the localization length £ seems to be 
promising, because it is sensitive to the nature, localized or 
extended, of the electron's eigenstate l34l l35tl . So far £ has 
been determined from the phase sensitivity of the ground-state 
energy II 1 71 12511 . Quite recently, Berkovits demonstrates that 
also the entanglement entropy can be used to extract the local- 
ization length ll36ll . However, in both methods the system size 
L should be always larger than £. 

Advantageously, the localization length can be extracted by 
a finite-size scaling analysis of the charge density structure 
factor even for L £, which works equally well for inter- 
acting systems l26ll . and therefore allows us to discuss the in- 
terplay between Anderson localization and CDW formation 
in a consistent manner. The charge structure is defined as 
0(q) - T2~:l=M n i+r) - (n{){n{ +r )]^\. Assuming 
an exponential decay of the equal-time density-density cor- 
relations in the Anderson insulating phase H^l, the structure 
factor scales with 



i 



2tt 2 



(4) 



e 8 « 



where q = 2ir/L <C 1 12611 . Equation (HJi contains two un- 
known parameters: the localization length £ and the disorder- 
modified TLL interaction coefficient K* . Hence, if the charge 
structure factor is determined numerically, £ and K* can be 
easily derived by fitting the numerical data with Eq. @. For 
vanishing disorder A — > 0, £ diverges and K* becomes the or- 
dinary TLL parameter K p . We are aware that a disordered ID 
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FIG. 2. (Color online) Charge structure factor of the AEM with A = 
0.1 and u) = 2 (left panel) and the AtVM with V/tf = 1.5 (right 
panel) sampled over 300 and 500 disorder realizations, respectively. 
Dashed lines give the finite-size scaling of C au {q) according to (0- 



system is no longer a TLL and consequently the TLL param- 
eter is ill-defined in the strict sense. Nevertheless, if the lo- 
calization length significantly exceeds the lattice constant, the 
short-range correlation functions should still show a power- 
law decay. Therefore we might gain some valuable informa- 
tion about the local motion of fermions from K*. 

Figure |2] demonstrates that the finite-size scaling of the av- 
eraged charge structure factor C a v(q) by means of (|4|i works 
to the best and equally well for the ID AEM and AtVM (this 
applies to all parameter values discussed below). To accomo- 
date the missing correlations owing to the OBC, we have plot- 
ted C av (<?) as a function of 1/ (L — 8) instead of 1/L (this way 
of plotting the data is nonessential but gives a quantitative re- 
finement of the fit). The parameter 6 is adjusted to reproduce 
K* = K p and £ = oo at A = 0. We note the general ten- 
dency that the charge correlations arising at finite L will be 
suppressed the more the larger A. 

In a next step, we extract the localization length £ and the 
modified TLL parameter for the disordered Edwards and t-V 
models. Figure [3] shows the dependence of £ and K* on the 
disorder strength A. First of all, we find a power law decay of 
£ with 1/A in the whole (A, luq; V) parameter regime: 




FIG. 3. (Color online) Left panels: Log-log plot of the localization 
length versus disorder strength for the AEM at fixed A = 0.1 (a), 
fixed u)o = 2.5 (c), and for the AtVM (e). Dashed lines are fits to 
Eq. with: (a) £ = 620, 98, 50 and 7 = 1.75, 1.65, 1.22 for 
A = 1, 0.2, 0.1; (c) £ = 500, 65, 40 and 7 = 1.81, 1.35, 1.1 for 
uo = 1, 2, 3; (e) £ = 440, 350, 230, 190, 150 and 7 = 1.95, 1.5, 
1.3, 1.1, 0.95 for V/tf = 0, 0.5, 1. 1.5, 2 (from top to bottom), 
respectively. Right panels: Corresponding results for the modified 
TLL parameter K* in the AEM [(b), (d)], and AtVM (f). 



£/£o = A -7 (5) 

[note that A is given in units of tb (</) in Fig.[3]for the AEM 
(AtVM)]. The estimated values of the (bare) decay length £0 
and the exponent 7 are given in the caption of Fig. |3]for char- 
acteristic model parameters. As expected, the localization 
length decreases with increasing disorder strength. Stronger 
electronic correlations, i.e, smaller A or larger ljq (larger V) in 
the AEM (AtVM), also tend to reduce £. In any case, £ turns 
out to be finite as soon as A > 0, indicating that the repulsive 
TLL, if realized for A = 0, makes way for an Anderson insu- 
lator. Thereby the localization length becomes comparable the 
lattice spacing at A = 2 in the AEM with A = 0.1, u> = 2.5, 
while it is still about 10 2 for the AtVM with V/t / = 2. 



The right-hand panels display striking differences in the A- 
dependence of K* for the models under consideration. These 
can be attributed to the fact, that the Edwards model con- 
tains two energy scales A and loq while the physics of the t-V 
model is merely governed by the ratio V/tf = 4/2 A, i.e,. 
uio drops out. Far away from the CDW instability, however, 
both models describe a weakly correlated TLL with K p < 1 
and K* slowly decays as the disorder A increases (see the red 
open triangles in Fig.|3j to make the comparison with the t-V 
model data easy, we have shown K* versus A/tf in the in- 
sets). If we move towards the CDW instability by decreasing 
A at fixed ujq > ujq, c or increase ujq at A < A c (cf. Fig. [U a 
non-monotonous behavior develops. At small A, K* is sig- 
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FIG. 4. (Color online) Left panels: Decay length £o as a function of 
the (effective) Coulomb repulsion Vt e f /) for the At VM and AEM. 
Right panel: Corresponding 7-exponent vs. K p , in comparison to the 
scaling relation l|6} (dashed line). For further explanation see text. 



nificantly enhanced as the disorder increases. Obviously weak 
disorder destabilizes the 2/cp-CDW correlations locally, since 
disorder-induced second- (and higher-) order boson-assisted 
(inelastic) hopping processes are possible in the AEM, even 
for wo 3> 1. This in sharp contrast to the AtVM, where only 
elastic scattering takes place and the intersite Coulomb repul- 
sion is hardly affected by A. As a result, in the disordered t-V 
model the CDW correlations will be stronger and more robust. 
Hence, for the AtVM, K* appears to be nearly independent 
from A for 0.5 < V/tf < 2. This also notably differs from 
the behavior found for the disordered Hubbard model l26ll . 
where the umklapp scattering is effectively enhanced by the 
formation of Mott plateaus appearing due to disorder 0811 . If 
A exceeds a certain value in the AEM, K* starts to decrease 
and finally the whole scaling procedure breaks down when 
£ > 1 (see the point at A = 2 in the upper right panel with K* 
well below 0.5). In this regime the wave functions of the par- 
ticles are strongly localized and the TLL-behavior is as much 
suppressed as the CDW correlations. Let us point out that the 
enhancement of K* triggered by the bosonic degrees of free- 
dom might serve as a explanation for the observed increasing 
charge velocity near a negatively charged defect in the single- 



wall carbon nanotubes 111 111 , since the TLL parameter K p is 



proportional to the charge velocity. 

We now focus on the localization behavior at large dis- 
tances [0(£ 3> 1)], and therefore made an attempt to an- 
alyze the decay length £0 an d the exponent 7, for both the 
AEM and AtVTVl, in terms of the interaction exponent K p 
and the charge susceptibility \c of standard TLL theory 13711 . 
We expect that £0 1S strongly influenced by the strength of the 
charge fluctuations quantified by Xc> which is given — for the 



t-V model — as \c 



2Kp 

7TV C 



arccos(- 



- 1 



with charge velocity v c . Figure [4] shows that £0 nicely scales 
with V/tf, i.e., £0 °c Xc m f act ( see upper panel). The same 
holds — perhaps surprisingly — for the AEM when £o/tj is 
plotted versus an effective intersite repulsion estimated from 
V e ff/tf = —2cos(ir/2K p ) (see lower panel). Moreover, for 
the disordered t-V model, the exponent 7 is connected to K p 
of the (spinless fermion) TLL system without disorder via the 



renormalization equation: d(A 2 )/dl = (3 
scale quantity I). This causes the scaling relation 1 37, 39, 4 



2K„)A 



(with 



7 = 2/(3 - 2K p ) . 



(6) 



The right panel of Fig.|4]displays that 7 basically depends on 
K p as predicted by Eq. ©. This means that the long-range 
localization properties of the AEM can be understood in the 
framework of AtVM with an effective intersite interaction in- 
duced by the bosonic degrees of freedom. Since the (effective) 
Coulomb repulsion tends in result to a lesser K p , 7 decreases 
with increasing Vt e ff\ (cf. Fig. |4j. In this way Anderson 
localization is suppressed by the 2fci?-CDW fluctuations trig- 
gered by V, which on their part also lead to charge carrier 
immobilization however. While 7 = 2 in the free-fermion 
limit (V, 1/A — > 0), it scales to unity approaching the CDW 
transition point, located, e.g., at A c ~ 0.07 for u>o = 2.5 re- 
spectively at wo, c — 3.1 for A = 0.1. 

The question how disorder affects the insulating CDW 
state could not be addressed by the above TLL-based scal- 
ing approach. Here, increasing the strength A of disorder, a 
crossover is expected from the CDW-regime dominated by the 
scale of the single-particle (charge) gap A c i(A = 0) 13 ill to 
the Anderson-regime, where A ^> A c i. In the former case the 
CDW survives moderate disorder (as experimentally observed 
for disordered Peierls-Mott insulators ilOll ). in the latter case 
a pseudogap scale may develop as in the strongly disordered 
Hubbard model |23fl. Figure [5] showing the spatial variation 
of the local fermion/boson densities for a specific but typi- 
cal disorder realization (note that any real experiment is per- 
formed on a particular sample), illustrates the situation deep 
inside the CDW phase (A -1 = 100, lo^ 1 = 0.4; cf. Fig.rjj. 
One realizes that long-range charge order ceases to exist but 
short-range CDW correlations may locally persist whenever 
neighboring onsite potentials will differ not much (see, e.g., 
the region i = 45 ... 55 in the lower panel of Fig. [5]). 

To summarize, using an unbiased numerical DMRG ap- 
proach, for the first time we investigated the interplay of disor- 
der and interaction effects including bosonic degrees of free- 
dom in the framework of the ID spinless fermion Anderson- 
Edwards model. Although the TLL phase disappears owing 
to the disorder, the localization properties of the Anderson in- 
sulator state can be understood in terms of scaling relations 
containing the charge susceptibility and the Luttinger liquid 
parameter of the metallic phase without disorder only, as in the 
case of the spinless fermion Anderson-t-V model. However, 
the Anderson-Edwards model reveals a more complex inter- 
relation between disorder and CDW correlations because ad- 
ditional scattering channels, involving bosonic excitation and 
annihilation processes, appear. This offers a promising route 
for adapting the description of low-dimensional transport in 
many disordered materials. Disorder also affects the CDW 
state in that true long-range order vanishes for strong on-site 
potential fluctuations though local CDW correlations survive. 

This work was supported by DFG through SFB 652. 
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FIG. 5. (Color online) Local densities of fermions (fjfi) (blue cir- 
cles) and bosons (&Jfo i ) (red open squares) in the central part of an 
open AEM chain without (A = 0) and with disorder (A = 2). Re- 
sults are given for a single realization {st} (black triangles). Model 
parameters are A = 0.01, lj = 2.5, and L = 128 (OBC). 
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